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ABSTRACT 

Context. Wavefront sensing in solar adaptive-optics is currently done with correlating Shack-Hartmann sensors, although the spatial- 
and temporal-resolutions of the phase measurements are then limited by the extremely fast computing required to correlate the sensor 
signals at the frequencies of daytime atmospheric-fluctuations. 

Aims. To avoid this limitation, a new wavefront- sensing technique is presented, that makes use of the solar brightness and is applicable 
to extended sources. 

Methods. The wavefront is sent through a modified Mach-Zehnder interferometer. A small, central part of the wavefront is used as 
reference and is made to interfere with the rest of the wavefront. 

Results. The contrast of two simultaneously measured interference- patterns provides a direct estimate of the wavefront phase, no 
additional computation being required. The proposed optical layout shows precise initial alignment to be the critical point in imple- 
menting the new wavefront- sensing scheme. 
Conclusions. 
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1. Introduction 

Adaptive-optic corrections for solar observations are currently 
based on Shack-Hartmann (SH) wavefront sensors, see e.g. 
Rimmele d2004b . Such sensors divide the wavefront by an array 
of lenslets, and each wavefront element forms a separate image 
on the detector. Local wavefront slopes are inferred from the im- 
age positions. When the source is a distant star, the images are 
small discs (either diffraction- or seeing-limited) and the image- 
centers are computed through a bary center calculation. For so- 
lar adaptive-optics, however, the image behind each lenslet is 
extended and the displacements must, therefore, be computed 
through a cross-correlation algorithm. 

With solar observations the number of lenslets is not lim- 
ited by flux considerations, as with the much fainter night-time 
targets, but by computational constraints. Even with the fastest 
available computers the calculations need to be optimized to 
meet the necessary correction-rates of ~ 2kHz. Thus, at the New 
Solar Telescope (NST) in Big Bear, California, the data from 
the 76 sub-aperture SH-sensor are analyzed with digital signal 
processors (DSP), programmed in assembly language, and the 
correlation is reduced to displacements within ±3 pixels, see 
e.g. Rimmele et al @04), Richards et al. (l2004l l2008t and 
Denker et al. (12007b . Suitable algorithms are currently being de- 
veloped for the planned upgrade to 308 sub-apertures. At the 
German Vacuum Tower Telescope (VTT) in Tenerife, Spain, the 
use of fast Fourier transforms permits correction frequencies of 
2100 Hz on a total of 36 sub-apertures. For the GREGOR tele- 
scope, Tenerife, Spain, the same team develops algorithms to 
sense the wavefronts at 2500 Hz on 156 sub-apertures (Berkefeld 
et al.. l20T0l) . 

A new wavefront- sensing technique is presented here. In this 
approach the phase is inferred without the need for a cross- 



correlation algorithm, and the computation requirements are 
greatly relaxed. A small, central part of the wavefront is used 
as reference and is made to interfere with the rest of the wave- 
front. A phase shift of n/2 is introduced in one arm of the inter- 
ferometer, and the two beams are then recombined via a beam- 
splitter cube. The wavefronts interfere with a phase difference 
from their common mean- value which varies around n/2 in one 
output, and -n/2 in the other output. The normalized intensity- 
difference between the two outputs is a direct measure of the 
wavefront phase. No additional computation is required. 

The method makes use of the Sun's brightness and is appli- 
cable to extended sources. Since all phases are measured relative 
to the same, central phase- value, the error propagation is inde- 
pendent of the number of sub-apertures. 

2. Formalism 

Let E be the electric field associated with the wavefront: 

E(x, y) = Eq • P(x, y) • exp(/ cp(x, y)) • exp(/ (kz - cot)) (1) 

The wave propagates along the z-axis in the orthogonal coordi- 
nate system (x,y,z). k = 2njX is the wave number, X the wave- 
length, co = 2nc/A denotes the angular frequency, c being the 
speed of light, cb is the wavefront phase and P the pupil trans- 
mission function. We consider a circular pupil of diameter D: 

D 2 

P(x,y) =1 if x 2 +y 2 < — 



- otherwise 



(2) 



In the proposed method, the wavefront is sent through a mod- 
ified Mach-Zehnder interferometer: the central part of the wave- 
front with diameter d = D/N is propagated along path (A), the 
complement is propagated along path (B). 
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- In path (A), an achromatic afocal system enlarges the central 
wavefront from diameter d to D - N d. Behind the afocal 
system the electric field equals: 



E A (x,y) = jfi-P(x,y) 

(x y \ 
i(p(—,—)j' exp (i (fa: - (ot)) (3) 

Path (B) introduces an achromatic n/2 phase- shift relative to 
path (A), and attenuates the amplitude of the electric field by 
a factor N 2 in order to make it comparable to path (A): 



E B (x,y) = jp-nx,y) 

• exp (i (0(jc, y) + n/2)) • exp (i (kz - cot)) (4) 

P' is the transmission function in path (B): P f (x,y) = 
P(x,y)-P(xN,yN). 

The two beams are recombined via a beam splitter, and the 
electric fields in the two outputs are then expressed by: 

E x (x,y) = —(E A (x,y) + E B (x,y)) 
V2 

E 2 (x,y) = -^-(E A (x,y) + exip(in)-E B (x,y)) (5) 
V2 

The additional phase-delay, n, for the field propagating from 
path (B) into output (2) is introduced upon reflection on the 
front of the beamsplitter mirror: the medium behind the mirror 
is glass, and has a higher refractive index than the air the field is 
traveling in. The field propagating from path (A) into output (1) 
reflects on the mirror-back, inside the glass of the beamsplitter 
and no additional phase shift is thus introduced. 
The resulting intensities are: 



E 2 

h(x,y) = j^'P\x,y) 
E 2 

I 2 (x,y) = jfi-P'(x,y) 



1 + 
1 + 



cos(0(x,>;) + 7r/2-0(^,^)) 

COs(0(x,j)-7T/2-0(^,^)) 



(6) 



These expressions do not apply inside the central disc of diam- 
eter d, where the field E B has zero amplitude and where both 
intensities equal: I\(x,y) = h(x,y) - E 2 /N 4 . In the following, 
only those points (x, y) are considered where P'(x,y) = 1. 

The central wavefront element - sent through path (A) - is 
used as reference. We therefore make the following approxima- 
tion: 



0(-,-)~0(0, 0) = 0, 
N N y 



So that, 



.,9 ? D 

if x +y < — 



(7) 



h(x,y) ~ j^.P'(x,y)[l-sin(cP(x,y))] 
E 2 

h(x,y) ~ j±.P'(x,y)[l + sm((f>(x,y))] 



(8) 



Under the condition of small phase distortions sin(0) ~ 0, 
and the intensity-contrast equals: 



h(x,y) - h(x,y) 
h(x,y) + I 2 (x,y) 



The phase at position (x,y) is inferred from only two in- 
tensity measurements, h(x,y) and h(x,y). The error propaga- 
tion from measurement to phase-estimate is thus independent of 
the number of measurement points (sub-apertures), while it in- 
creases logarithmically with the number of sub-apertures when 
SH-sensors are employed (Hudgin |l9771 Fried l 19771 Kellerer & 
Kellerer[20ni). 

Turbulence-induced path-length distortions, S, are achro- 
matic, hence phase distortions are inversely proportional to the 
wavelength: cp(x,y, A) = 2nd(x,y)/A. If the spectral range of the 
phase sensor is [Aq, Aq + AA], Eq.[9]becomes: 



h(x,y) - h(x,y) 2n 
h(x,y) + I 2 (x,y) 



— .ln(l + AA/Ao).6(x,y) 

A/t 



(10) 



(9) 



Wavefront sensing in solar adaptive-optics is typically done in- 
side a 50 nm wide spectral-band centered around 550 nm, in this 
case Z i t^(-l) z (A/i//io)7(*+ 1) ~ 1. This band is assumed for our 
sensor. 



3. Implementation of the phase sensor 

3.1. Optical Layout 

Fig. Q] sketches the optical layout as it would be used in the case 
of a point source. Along path (A) the central part of the wave- 
front is enlarged from its initial diameter d to the pupil diameter 
D = Nd. Along path (B) an achromatic phase-retarder intro- 
duces a 7i / 2 phase-shift relative to path (A). An absorption plate 
attenuates the field-amplitude by a factor N 2 in order to approxi- 
mately equalize the amplitudes in the two paths. Two delay lines 
compensate for the increment in path-length introduced by the 
achromatic afocal system in path (A). 

The beams are recombined via a beam-splitter. For a flat in- 
coming wavefront, the phase difference upon recombination is 
7i 1 2 in output (1) and -n/2 in output (2). For a distorted wave- 
front - and if the phase variance of the reference wavefront 
(path A) is negligible - the phase difference upon recombina- 
tion equals (cp(x,y) + n/2) in output (1) and (cp(x, y) - n/2) in 
output (2). 

The mean-phase difference introduced by a path-length dif- 
ference, dL = L B - L A , between paths A and B equals: dcp =< 
cp B > - < cpA >= 2n - dL/A. A path-length difference controlled 
within ~ A/20 ensures that the instrumentally induced phase dif- 
ferences are well below n/2, and that the phase differences vary 
around n/2 and -n/2 in the two outputs. 

The method is not actually intended for application to point 
sources because, on the one hand, the intensity attenuation factor 
1 /N 2 would then not be acceptable and, on the other hand, be- 
cause any remaining phase fluctuations within the central wave- 
front segment may, with point sources, make the broadened 
beam in path (A) insufficiently uniform. As will be seen, both 
limitations are less critical in solar observations. 



3.2. Adjustment for an extended object 

To allow for the image extension in solar observations, the op- 
tical layout is adjusted so that the interference patterns from 
sources in different directions will be properly superimposed 
(see Fig.[2|: 
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from telescope p^thi (A) 




Fig. 1. Optical layout of the phase sensor. Path (A) propagates the central part of the wavefront. The achromatic afocal system 
(Li, fl , Li, a ) enlarges the wavefront diameter from d to D = N d. The two delay-lines compensate for the path-length increment due 
to the afocal system in path (A). 



- Mirror M\ , which sends the central part of the wavefront into 
path (A), is placed in a pupil plane, where the wavefronts 
from different angular-directions overlap. This ensures that 
the central parts of all wavefronts are picked up by M\ . 

- The afocal system (L\ A ,Li,a) is adjusted to re-image the 
pupil plane (Mi) onto the two detectors (P\ and P2). This 
requires L\ A and Li, a to be convergent lenses, and M\ to 
be placed between L\ A and the focal plane of L\ A . If L\ A 
were divergent and La, a convergent, the image of M\ through 

Li^a) would be virtual, and could, thus, not be projected 
onto the detector planes, P\ and P2. 

- Another afocal system (L^, Lzj) is placed in path (B) to re- 
image the pupil-plane, M\, onto the two detectors, P\ and 
Pi. 

With the adjusted optical layout, the measured contrast 
equals the intensity weighted phase-average over a solid-angle, 
Q: 

Idx,y)-I 2 (x 9 y) ff a m<Kx 9 y 9 0)d 2 O 
h (x, y) + I 2 (x, y)~ fj a L(0) d 2 6 

L(6) [m~ 2 s" 1 sr" 1 ] is the solar radiance from the angular direc- 
tion e. 

Adaptive optics with a single deformable mirror is in princi- 
ple limited to acceptance angles within the isoplanatic diameter 
(a = 10" for typical daytime observations). For larger accep- 
tance angles, Q > na 2 and the correction can merely account 
for the mean phase shift over Q. Since this mean value is pre- 
dominantly determined by the nearby turbulence, the resulting 
approximate correction is termed ground-layer adaptive optics. 
A second condition for exact adaptive optics is the sufficiently 
close spacing of actuators. Their separation should correspond 
to not more than the Fried distance (roughly 0.05 m for typical 



day-time observations), otherwise the correction accounts again 
merely for a mean phase shift. 

With a point object, the central sub-aperture contains, de- 
pending on its size, a certain level of phase differences. Ideally 
one would wish to average out these differences in the widened 
beam (A) to make it as uniform as possible. While this can not 
be done, a somewhat increased angle of acceptance will have the 
similar effect of replacing the phase shifts at each point by a lo- 
cal average which makes the central beam more uniform. The 
proposed method is, thus, inherently suited for extended images. 

4. Sensitivity of the phase sensor 

4.1. Maximum apparent magnitude 

Assume that sensor areas equal the size of the central sub- 
aperture. The broadened beam from the central sub-aperture is 
then required to interfere with the N 2 - I other areas. The num- 
ber, Z, of photons per detector read-out must, thus, be sufficiently 
high: 

Z = I S At>Z 'N 2 (12) 

/ [m~ 2 s" 1 ] is the target's irradiance, S the collecting surface of 
the central sub-aperture and At is the exposure time. Zo denotes 
the minimum number of photons for the interference in each sub- 
area. For Zo = 100 the signal to photon-noise ratio stays below 
10%. 

Initial calculations (see Section |4~2b suggest that the root- 
mean-square (rms) phase-deviation, cr^, over the whole wave- 
front needs to remain within ~ 1.2 rad. The central part of the 
wavefront can serve as a reference if its rms phase-deviation, cr , 
is substantially smaller than cr^, say cr < cr^/5. Accordingly, 
the diameter of the region on the telescope aperture that corre- 
sponds to the central sub-aperture should ideally not exceed a 
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Fig. 2. Adjustment of the optical layout (Fig.[T]) for an extended source such as the Sun. Both afocal systems, (Li^ a ,L2, a ) and 
Li,b), re-image the pupil plane, Mi, onto the detector planes, P\ and P^- 



fraction of the Fried length, ro (Roddier 1981): 
^<ro-<i x = 0.19.r 



(13) 



intensity difference ceases to be proportional to the phase distor- 
tion. Instead (see Eq.[2]): 



With the diameter of the NST telescope, D = 1.6 m, and a Fried 
length r = 0.05 m, we then get: d - 0.01m and N = D/d = 
168. Daytime adaptive-optics correction is typically done at a 
frequency of / = 1 /At = 2 kHz, so that: 

/> 8- 10 13 m- 2 s _1 (14) 
And the target's apparent-magnitude needs to be below: 

m < -2.5 log - = -10.5 (15) 

Iq = 5 • 10 9 m~ 2 s" 1 equals the irradiance from a star with zero 
apparent magnitude inside a 50 nm wide spectral-band centered 
around 550 nm (Allen's Astrophysical Quantities, 2000). 

When the wavefront sensing is done inside the isoplanatic 
cone (typically a = 10"), the corresponding apparent magnitude 
of the Sun equals: m = -26.7 + 2.51og(0.5 • 3600/10) = -21.1. 
The number of photons received in the central sub-aperture is 
thus by far sufficient. 

4.2. Maximum amplitude of the phase-distortions 

In exploratory computations different 288 x 288 matrices of 
phase values have been generated with Kolmogorov statistics, 
and the normalized intensity-differences, (I\ - h)l(h + h), have 
been calculated in terms of Eqs.[3]-[9] The central 32 x 32 el- 
ements were taken to represent the reference wavefront that is 
propagated through path (A) of the interferometer. The normal- 
ized intensity differences are shown in Figs.[3]-|5] for different 
phase screens, with rms phase-deviations between cr^ = 0.1 rad 
and 1.2 rad. cr^ = 1.2 rad appears to be a limit: for larger phase 
distortions, sin(^) ~ fails to be a valid approximation, and the 



h(x,y) - I 2 (x,y) 
7i(x,v) + / 2 (x,v) 



sin(0(x,y)) 



(16) 



The phase values need then to be unwrapped, and the sens- 
ing scheme looses its advantage of minimal computational re- 
quirements. For a D = 1.6 m telescope and a Fried length 
ro = 0.05 m, the turbulence-induced rms phase-distortions equal 
in fact: cr = (D/r ) 5/6 = 18 rad (Roddier|l98B. One must, 
therefore, take into account that the phase differences between 
beams (A) and (B) are determined only modulo 2n which re- 
quires an initial computational procedure which might start out 
the adaptive-optics correction-loop on a central subset of the sen- 
sor field and progressively extend the set as the correction takes 
effect, this extension being either algorithmically coded or be- 
ing controlled by operating a diaphragm in the collimated beam. 
Once an approximate phase correction is achieved the precise 
correction can proceed without further computational effort at 
the desired frequency. 



5. Conclusion 

A scheme is presented to sense wavefront distortions on bright, 
extended sources. The central part of the wavefront is used as 
reference and is made to interfere with the rest of the wavefront 
by use of a suitably modified Mach-Zehnder interferometer. The 
intensity contrast in the two outputs of the interferometer is a 
direct measure of the wavefront-phases. Once an approximate 
phase adjustment is attained the continued exact correction re- 
quires no additional calculation, i.e. the computation require- 
ments are greatly relaxed in comparison to the current use of 
correlating SH-sensors. The most critical step in the implemen- 
tation of the new wavefront- sensing method is likely to be the 
initial alignment of the Mach-Zehnder interferometer, since the 
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Fig. 4. Same as on Fig.O with a different realization of phase values. The rms phase deviation equals 0.6 rad. 
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Fig. 5. Same as on Fig.[3j with a different realization of phase values. The rms phase deviation equals 1.2 rad. 



path lengths in both arms of the interferometers need to be equal- 
ized within A/20 (see Section l3Jl) . This, as well as the suitable 
control strategy will need to be explored in an experimental set- 
up. 

The phase differences are deduced from the contrast of two 
simultaneously measured intensities and, accordingly, the esti- 
mates are not affected by scintillation. Furthermore, all phases 
are measured relative to the same, central phase- value, and the 
error-propagation factor is therefore independent of the number 
of sampling elements. This contrasts with the logarithmic in- 
crease of the propagation factor when increasing numbers of SH 
sub-apertures are used. 

The method requires the flux in the central sub-aperture to be 
sufficiently high, and is thus especially suited for solar adaptive- 
optics. 
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